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Plasma ionized calcium (Ca ) concentrations are tightly regulated in the body and maintained within 
a narrow range; thus it is challenging to quantify calcium absorption under normal physiologic 
conditions. This study aimed to develop a mechanistic model for the parathyroid hormone (PTH) 
response after calcium intake and indirectly compare the difference in oral calcium absorption from 
PTH responses. PTH and Ca 2+ concentrations were collected from 24 subjects from a clinical trial 
performed to evaluate the safety and calcium absorption of Geumjin Thermal Water in comparison 
with calcium carbonate tablets in healthy subjects. Indirect response models (NONMEM Ver. 7.2.0) 
were fitted to observed Ca 2+ and PTH data, respectively, in a manner that absorbed but unobserved 
Ca 2+ inhibits the secretion of PTH. Without notable changes in Ca 2+ levels, PTH responses were 
modeled and used as a marker for the extent of calcium absorption. 

Key Words: Calcium, Indirect response model, Nonlinear mixed effects modeling, Parathyroid hormone 



INTRODUCTION 

Calcium levels are tightly regulated in the body. In ad- 
equate calcium homeostasis, ionized calcium concentrations 
([Ca 2+ ]) are maintained within a narrow range (1.1 ~ 1.35 
mmol/L) [1], and any derangement from it causes serious 
health problems. Parathyroid hormone (PTH), synthesized 
and secreted by the chief cells of the parathyroid glands, 
is one of the key endocrine regulators of extracellular cal- 
cium concentrations [2,3]. Calcium-sensing receptors, lo- 
cated on parathyroid cells, are activated by either increases 
or decreases in [Ca 2+ ], which result in the inhibition or stim- 
ulation of PTH secretion, respectively [2,3]. In turn, serum 
PTH acts on bone formation/resorption, and Ca 2+ absorption 
and excretion from the bone, kidney, and intestine to main- 
tain [Ca 24 ] at the normal physiological range [2,3]. 

In order to evaluate the safety and calcium absorption 
of Geumjin Thermal Water, a randomized, parallel clinical 
study was conducted. Due to Ca 2+ homeostasis, it is chal- 
lenging to quantify oral calcium absorption unless a more 
sophisticated approach is used, such as tracer methods (i.e., 



Received January 30, 2014, Revised February 21, 2014, 
Accepted March 14, 2014 

Corresponding to: Dong-Seok Yim, Department of Pharmacology, 
College of Medicine, The Catholic University of Korea, 222 Banpo- 
daero, Seocho-gu, Seoul 137-701, Korea. (Tel) 82-2-2258-7327, (Fax) 
82-2-2258-7876, (E-mail) yimds@catholic.ac.kr 




This is an Open Access article distributed under the terms of the 
Creative Commons Attribution Non- Commercial License (http:// 



creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial 
use, distribution, and reproduction in any medium, provided the original work 
is properly cited. 



the use of a suitable calcium isotope). However, such meth- 
ods are not readily applicable to most marketed calcium 
products. 

In this study, it was hypothesized that PTH levels would 
carry more information about calcium absorption than the 
Ca 2+ concentration itself, as any increment in plasma [Ca 2+ ] 
after the intake of calcium would have resulted in a change 
in PTH levels, which is not prone to be tightly regulated. 
In order to utilize all of the information included in the 
PTH responses that are changed by time, more sophisti- 
cated, model-based approaches are preferred to a simple 
statistical test. K-PD (K for kinetic instead of PK for phar- 
macokinetic, and PD for pharmacodynamic) modeling [4] 
is currently used to model PD response-time profiles when 
the drug concentrations are not available. Similar to the 
rationale for K-PD modeling, the PTH response was also 
modeled without noticeable changes in observed calcium 
levels. Specifically, the goals of this analysis were to devel- 
op a mechanistic model for PTH and to compare the degree 
of calcium absorption after ingestion of thermal spring wa- 
ter or a calcium carbonate supplement using the PTH 
response. 

METHODS 

Study design 

This was a randomized, parallel clinical trial of Geumjin 
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Thermal Water (thermal spring water hereinafter) to eval- 
uate its safety and calcium absorption in healthy Korean 
subjects. The safety and the extent of calcium absorption 
of thermal spring water were compared with that of a cal- 
cium carbonate supplement. The inclusion criteria in- 
cluded: 1) healthy adults between 20 and 50 years of age 
at screening, 2) male and female subjects whose weights 
were greater than or equal to 55 kg and 45 kg, respectively, 
and within±20% of ideal body weight, calculated as (height 
in cm — 100)x0.9, and 3) subjects must be willing and able 
to comply with scheduled visits, treatment plan, laboratory 
tests, and other trial procedures after being informed about 
the study and agree to provide a written consent. Subjects 
allergic to calcium supplements, whose 25(OH) Vitamin D3 
levels were outside the range of 4.8 ~ 52.8 ng/mL at screen- 
ing, whose calcium levels were outside 8.0 ~ 10.0 mg/dL, or 
otherwise were considered unhealthy according to the study 
protocol were excluded from the study. The demographic 
summary of the study population is provided in Table 1. 

Twenty-four subjects were randomized to one of three 
treatments (2:1:1 ratio). During the mornings of days 
1 and 7, 12 subjects received 240 mL of Geumjin thermal 
spring water that contains 400 mg of elemental calcium, 
followed by 100 mL of purified water; 6 subjects received 
500 mg of calcium carbonate (CaCOs) (two tablets contain- 
ing 200 mg calcium each) with 240 mL of normal saline 
plus 100 mL of purified water; and 6 subjects received 500 
mg CaC03 with 340 mL of purified water. On days 2 
through 6, the morning (daily) doses described above were 
split into two doses and administered twice daily. Blood 
samples were collected before and up to 8 hours (0, 0.5, 
1, 1.5, 2, 3, 4, 6, and 8 hours) after the calcium intake. 

Model development 

While the PTH-Ca 2+ system is very complex, it was sim- 
plified with the following assumptions: 1) the net absorbed 
calcium from each calcium supplement is positive but not 
observed due to the Ca 2+ homeostasis; 2) the drop in the 
PTH response is only due to the increased but unobserved 
Ca 2+ , as a negative feedback; and 3) no circadian rhythm 
nor systemic differences in baseline PTH levels on days 1 
and 7 were present, other than random variability. The 
schematic of the model is presented in Fig. 1. Indirect re- 
sponse models with zero-order input and first-order output 



rates were fitted to the observed Ca 
spectively, as follows. 
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Table 1. Subject Demographics 



Mean (range) 



Age (years) 
Weight (kg) 
Height (cm) 
Sex (Male/Female) 
Smoking (Yes/No) 
Alcohol (Yes/No) 
Caffeine (Yes/No) 



26 (21-39) 
68.5 (55.1-79.3) 
174.5 (161-184) 
22/2 

7/17 
14/10 
14/10 
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Where kinc a represents the zero-order input rate for Ca , 
koutca represents the first-order output rate for Ca 2+ , 
kinpiH represents the zero-order input rate for PTH, and 
koutpiH represents the first-order output rate for PTH. 

Increased Ca 2+ levels due to the intake of thermal spring 
water or the calcium supplement were believed to decrease 
the secretion of PTH, which was incorporated into the mod- 
el as an inhibitory Emax model as follows. 
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Where Chi 



is the net absorbed but unobserved Ca 



ka represents the first-order absorption rate constant, 
Emax represents the maximum effect of [ Ch u * gbs ], on PTH 
secretion (fixed to 1), 

EC50 represents the concentration of [Ch 2 u + obs ] to achieve 
half of Emax, while other terms are as defined previously. 

[Ctfiepot] a< t=0 was assumed to be 1. 

The bioavailability of calcium carbonate was estimated 
relative to that of thermal spring water. The type I error 
rate (a) for covariate testing was set as 0.05. Since the 
number of subjects was small and the main interest was 
the difference between thermal spring water compared with 
calcium carbonate, the 12 subjects that received calcium 
carbonate tablets were combined and treated as one treat- 
ment (i.e., calcium carbonate). No other covariates than the 
formulation of calcium were tested. Model development was 
guided by the objective function value, goodness of fit, and 





Fig. 1. Schematic of the model, ka, first-order absorption rate 
constant for Ca; F, relative bioavailability; kin_ca, zero-order input 
rate for Ca 2+ ; kout_ca, first-order output rate constant for Ca 2+ ; 
kin_pth, zero-order input rate for PTH; kout_pth, first-order output 
rate constant for PTH. *Unobserved. 
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the physiological plausibility of the model and parameter 
estimates. 

Model evaluation 

A basic set of diagnostic plots was generated for both Ca 2+ 
and PTH. For a visual predictive check, 500 data sets with 
the same study design of the observed data set were simu- 
lated using the final model and parameter estimates, and 
the 5-50-95 th percentiles were compared with those ob- 
served in this study. For bootstrap, 500 data sets were sam- 
pled with replacement and nonparametric 95% confidence 
intervals (2.5 ~ 97.5 percentiles), and the median values 
were assessed using the runs with successful minimization. 



Software/ platform 

The analysis was performed using NONMEM 7.2.0 
(ICON Development Solutions, MD) in the Pirana 2.6.1 
platform [5]. R 2.15.1 (www.r-project.org) with R studio ver- 
sion 0.96.330 (http://www.rstudio.com/) was used for data 
formatting and plotting. Perl-Speaks-NONMEM 3.5.3 was 
used for bootstrapping. The estimation method used was 
the first-order conditional (FOCE) with interaction, while 
ADVAN13 was used to solve differential equations. 



RESULTS 

A total of 423 PTH and 423 Ca 2+ concentrations, collected 
from 24 subjects, were included in the analysis. The final 



Table 2. Parameter Estimates, 95 % Confidence Intervals (CI), and Bootstrap Results 



Parameters 


Estimate 


95% CI 


Bootstrap median 


95 lh Percentiles 


kin_ca 


3.41 


3.39, 3.43 


3.41 


3.39, 3.43 


kout_ca 


2.86 fixed 








kin_pth 


21.6 


12.4, 30.8 


20.6 


9.12, 33.0 


kout_pth 


0.849 


0.513, 1.19 


0.821 


0.375, 1.26 


ka (hr-1) 


0.796 


0.438, 1.15 


0.871 


0.54, 7.18 


EC50 (mmol/L) 


0.158 


0.0924, 0.224 


0.148 


0.026, 0.232 


Relative Fl 


1.98 


1.06, 2.90 


2.14 


1.31, 8.58 


IIV in kin_ca 


0.000229 


0.000814, 0.000377 


0.00021 


0.0001, 0.0004 


IIV in kin_pth 


0.0453 


0.0260, 0.0646 


0.0388 


0.0213, 0.0616 


IIV in ka 


0.406 


0.161, 0.651 


0.390 


0.164, 1.34 


IOV in Kin_pth 


0.00969 


0.0015, 0.0179 


0.009 


0.0026, 0.0189 


SD_ca (mmol/L) 


0.0258 


0.0212, 0.0304 


0.0259 


0.022, 0.0301 


CV_pth (% CV) 


21.0 


18.8, 23.2 


20.9 


18.9, 23 



IIV, inter-individual variability; IOV, inter-occasion variability; SD, standard deviation; CV, coefficient of variation. 



Ionized Calcium 




0 2 4 6 3 1.16 1.18 1.20 1 22 



Fig. 2. Diagnostic plots for observed Ca + . DV, observed Ca + ; PRED, model predicted values; IPRED, individual model predicted values; 
CWRES, conditional weighted residuals. Gray lines are loess smooth. 
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parameter estimates and 95% confidence intervals are pro- 
vided in Table 2. The elimination rate constant for Ca 2+ 
was fixed to the literature value [6] of 2.862 hr , as both 
the calcium turnover rate and the fractional rate of elimi- 
nation might not be identifiable unless some extreme 
change (such as chelation of extracellular calcium by EGTA, 
ethylene glycol tetraacetic acid) is introduced into the 
system. The same constant was used for the rate of elimi- 
nation for unobserved calcium in the model (Fig. 1). The 
first-order elimination rate constant for PTH was estimated 
to be 0.849 hr 1 , which was smaller than the values found 
in the literature (as a half-life, less than several minutes 
[6,7]). The first-order rate constant for calcium absorption 
was estimated to be 0.796 hr , which could have been in- 
fluenced by the rate of calcium absorption as well as the 
change in the PTH response in this type of model. It was 
also necessary to fix the Emax to 1, as the maximal calcium 
effect on PTH was unlikely to be observed with the calcium 
amount evaluated in this study; yet it was assumed that 
100% suppression in PTH secretion can be achieved if the 
maximal calcium dose is introduced. 

In order to account for the difference in baseline PTH 
level, inter-occasion variability (IOV) was included into the 
zero-order input rate constant for PTH. In addition to the 
random individual difference of 21.4%, 9.84% IOV was esti- 
mated in this study (Table 2). 

The relative bioavailability of calcium carbonate was esti- 
mated to be 1.98 (RSE 24%) times higher than that of ther- 
mal spring water (p<0.05). The treatment effects on the 
calcium absorption rate constant (ka) or on EC50 were not 
significant (a =0.05). 

Basic diagnostic plots, such as observed (DV) versus 
(individual) predicted values (I/PRED), showed that the fi- 
nal model described well both the observed Ca 2+ and PTH 
levels (Figs. 2 and 3, respectively, upper panel). As sup- 
ported by the even scattering in the conditional weighted 



residual (CWRES) plots versus time or individual predicted 
values (Figs. 2 and 3, lower panel), no apparent bias or sys- 
tematic trend seemed to exist. 

A visual predictive check revealed the congruency be- 
tween the observed data and the model predictions at each 
occasion (i.e., at Days 1 and 7) (Figs. 4 and 5, upper panel) 
or per treatment (Figs. 4 and 5, lower panel) for both Ca 2+ 
and PTH. 



DISCUSSION 

From a clinical trial to evaluate the safety and calcium 
absorption of Geumjin thermal spring water in comparison 
with calcium carbonate tablets, the PTH and Ca 2+ levels 
were measured for 8 hours after oral administration. Daily 
doses of thermal spring water or calcium carbonate tablets 
(either with normal saline or drinking water) contained 400 
mg of elemental calcium. When the raw data were visually 
inspected (not shown), it was clear that the PTH response 
was a more sensible indicator than Ca 2+ itself; therefore, 
the PTH response was used to evaluate calcium absorption 
in this analysis with the assumptions described in the 
Method section. The rationale was that the net absorbed 
calcium from oral intake could have inhibited the secretion 
of PTH, but the calcium itself is not noticeable due to cal- 
cium homeostasis. If the bioavailability of calcium was sim- 
ilar between thermal spring water and calcium carbonate 
tablets, it was assumed that the change in PTH would have 
also been similar. When the treatments - thermal spring 
water versus calcium carbonate tablet - were tested for rel- 
ative bioavailability, it appeared that calcium carbonate 
had an approximately two times higher bioavailability than 
thermal spring water. However, the finding needs to be 
carefully interpreted, as the test was performed on hypo- 
thetical calcium levels, not the actual observations. Because 
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Fig. 3. Diagnostic plots for observed PTH. DV, observed Ca + ; PRED, model predicted values; IPRED, individual model predicted values; 
CWRES, conditional weighted residuals. Gray lines are loess smooth. 
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Fig. 4. Visual predictive check for Ca + . Obs, observed data; Sim, simulated values. 



baseline PTH changes by circadian rhythm were not in- 
corporated into the model, the exact pattern of PTH change 
by calcium absorption is not known in this study. 

It is worthwhile to mention previous efforts to model the 
PTH-Ca 2+ system. Abraham et al. [6] also used indirect re- 
sponse models, but they were able to distinguish the syn- 
thesis of preproparathyroid hormone and the secretion of 
PTH by integrating human data and rat experiments where 
EGTA was infused to chelate plasma Ca 2+ . Shrestha et al. 
[8] developed a mathematical model for the PTH response 
to acute changes in Ca 2+ . They characterized the short-term 
dynamic of PTH response and Ca 2+ concentrations as the 



reverse sigmoid relation and estimated parameters using 
clinical data collected from three subjects who underwent 
hypocalcemic clamp tests [8]. On the other hand, our study 
focused on estimating the relative degree of calcium absorp- 
tion rather than developing sophisticated feedback models. 
Although the mild condition in our study would limit the 
implementation of more complex and mechanistic models, 
we believe a step was taken forward by fully utilizing the 
time course of PTH responses, rather than only the change 
from the baseline. Heaney [9] calculated the area under the 
curve (AUC) of Ca 2+ over 9 h by measuring radioactivity 
after giving radiolabelled calcium carbonate formulations. 
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Fig. 5. Visual predictive check for PTH. Obs, observed data; Sim, 



His report is unique in that the change in the serum Ca 
concentration was directly measured; however, radiolabe- 
ling is not easily applied to human studies due to the cost 
and time involved. 

In this study, semi-mechanistic modeling was performed 
to characterize the Ca 2+ -PTH system after intake of thermal 
spring water containing calcium or calcium carbonate 
tablets. Without noticeable differences in the plasma Ca 2+ 
levels, the changes in the PTH responses were used as a 
surrogate marker for calcium absorption. With a more en- 
hanced study design, it is believed that our approach can 
be applied to evaluate calcium absorption, as determined 
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simulated values. 



by PTH-time responses. 
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